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1. Introduction 

Response of a many-body system to external fields provides various information about 
intrinsic dynamics of the system. One of the efficient approaches to get a deep 
insight into a nuclear structure is to study a response of a nucleus to an external 
Coriolis field. In particular, we can trace the interplay between single-particle and 
collective degrees of freedom by studying a rotational dependence of the kinematical 
and dynamical moments of inertia. In its turn, these quantities are benchmarks for 
nuclear microscopic models, that allow to understand main features of a nuclear field. 
For example, the importance of pairing correlations, introduced in nuclear physics by 
Bogoliubov £Q and Bohr, Mottelson and Pines :2. in an analogy with the correlations 
in superconductors, had been recognised first in the description of nuclear inertial 
properties [Hj- 

The kinematical moment of inertia is described quite well within the cranking 
approach, for instance, in the relativistic Hartree-Bogoliubov and non-relativistic 
density-dependent Hartree-Fock-Bogoliubov theory with Gogny forces (see 0j and 
references therein). However, there is a systematic discrepancies between those 
theories and experiment with regard of the dynamical moment of inertia. It is quite 
desirable to understand the main source of discrepancies between experiment and 
theory, since the above approaches represent state of art of nuclear structure studies. 

It is well understood that the two-body correlations are important for a correct 
description of the moments of inertia (see discussions in 4 ) . One of the systematic 
ways to include these correlations is the random phase approximation (RPA) based 
on a self-consistent solution of the cranking mean field equations |5l E] ■ This enables 
one to also restore symmetries broken on the mean field level (see textbooks pOIH])- 
Practical realisation of these aspects has been done recently in a cranking harmonic 
oscillator model with a residual quadrupole-quadrupole interaction [5]. In particular, 
it was shown that the Thouless-Valatin moment of inertia JU| (which, in general, 
contains Belyaev [Hj , Migdal and other terms) calculated in the RPA is equivalent 
to the dynamical moment of inertia calculated in the mean field approximation 



Here 51 denotes the rotational frequency, is a self-consistent solution of the mean 
field equations and Emf is a total mean field energy in the rotating frame. 

Quantum oscillations around the mean field solutions provide the additional 
contribution to the total energy. How such a quantal effect can be obtained within the 
RPA for the case of rotating systems is generally described in |12l IT3] . However, the 
practical application of the RPA in rotating deformed nuclei is a demanding numerical 
task. Therefore, the influence of quantum fluctuations on the moment of inertia is 
till now scarcely studied. With regard of the easier case of pairing fluctuations more 
extensive investigations were performed (c.f. The first attempt to incorporate 

the pairing and quadrupole oscillations |I6| was done in a restricted configuration space 
and the self-consistency between the residual interaction and the mean field was not 
taken into account. A self-consistent analysis carried out in a simplified model [5] 
demonstrated that the quantum correlations modify the Thouless-Valatin moment of 
inertia. Therefore, there is a strong motivation to understand the role of quantum 
correlations in realistic calculations for inertial properties. 

The aim of the present paper is to study exclusively the influence of shape 
vibrations on the angular momentum and on the moment of inertia. Starting from the 



'Emf 
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self-consistent cranking mean field calculation the total binding energy as a function 
of the rotational frequency has to be calculated up to RPA order, i.e. with inclusion 
of the RPA correlation energy. Our calculations are based on separable QQ-forces as 
an effective interaction. The mean field part of such calculations could be carried out 
with more advanced effective interactions like e.g. Skyrme or Gogny type forces. 

However, the RPA correlation energy is not feasible for non-separable forces 
because of the too large configuration space needed for the cranking model. In 
addition, the problem of spurious solutions in the RPA calculations for the state of 
art effective interactions is silently avoided in the literature even for a zero rotational 
frequency. In other words, again the question arises about a consistency between 
a mean field calculations and the RPA contribution for those type of forces. On 
the other hand, the QQ-forces are the most important interaction term for inducing 
nuclear deformation and rotational excitations. The mean field plus RPA calculations 
with the QQ-forces based an a realistic cranking mean field potential with a spin-orbit 
coupling are by no means trivial. As we shall describe (see below) such calculations 
can be executed in different ways which enable us study specific properties of these 
forces and their effects in more detail. 

We will consider the following versions of the QQ-model: 

(i) the QQ-forces with a full N-mixing of oscillator shells (AN = 0, ±2); 

(ii) the QQ-forces without N-mixing (AN = 0) which is the commonly used a QQ- 
model (e.g. |7]); 

(iii) thc double-stretched QQ-forces with the volume conservation constraint |17l I18j . 
By comparing results of the various treatments (i-iii) provides new insights in the 
nature of both the QQ-forces and their influence onto the rotational properties. 
Because of the schematic character of the QQ-forces our investigations are merely an 
exploratory approach to the influence of quantum shape fluctuations upon rotational 
properties rather than a quantitative description of experimental data. The latter 
would require more extensive numerical efforts that are not intended in this paper. 

The paper is organised as follows. In Sec. we introduce our mean field and RPA 
models. The influence of the volume conservation constraint on the RPA correlations 
is discussed in Sec. EH This is followed by a summary and discussion in Sec. 0] 

2. The Model 

Our mean field plus RPA calculations are based on the following Hamiltonian: 

^ = E efe 4^-f E Q 2 m = H sph + H QQ . (2) 

k m=0,±l,±2 

Here are the single-particle energies of the spherical oscillator Hamiltonian h sp h 
? p 2 M ,„ 2 

fi.ni, = 1 ujr.r . (3 

p 2M 2 w 

The operators c\. (£k) are creation (annihilation) operators with the suffix fc labeling a 
set of quantum numbers. For the sake of convenience we have chosen in equation 
a quadratic form of the QQ-forces by taking a set of Hermitean quadrupole operators 
Q m built up from r 2 Y2o and the linear combinations f 2 (Y2 m ± ^2-m) for rn = (1,2). 
The operators Q m read as 

Qo = V^o E 9o,kl (4^ + 4a) . ( 4 ) 

kl 
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with 



M {c\ci±c\c^j , (5) 
ci ± ctc r ) (6) 



, 1 to = 0, —1, 2 
<*"H 1 TO=l-2 ^ 




where the index k{l) is labeling a complete set of quantum numbers in order to form 
the matrix elements q m ,M = (k\Qm\l)- 

The index k refers to the time conjugated state. The sums in equations 
split in a proton and neutron part. We consider the isoscalar quadrupole interaction 
only, since isovector terms can be treated in a similar way. Hereafter, any isospin 
index is omitted. 

The Hamiltonian H , equation J5J is, of course, formally the same expression for 
all the versions (i-iii) of the QQ- model described in the introduction. What is different 
concerns the particular treatment of the Q-operator (N-mixing or not and stretching 
transformation), the adjustment of the spherical oscillator frequency, u>o, to account 
for the volume conservation condition and the adaption of the strength parameter k. 
This is because the rotational invariance of the total Hamiltonian is required for a 
clean treatment of the angular momentum conservation within the RPA |12l ll.'ij . 

We consider the Routhian 

H' = H - VLJ Z (8) 

to describe the rotational properties of the system. Here J z is the angular momentum 
component about the z-axis. Note that the cranking term in H' does not change the 
exact wave functions but only fixes the angular momentum about the quantisation 
axis (z). 

The relevant mean field part of the Routhian H' is originated from the 
Hamiltonian H in equation J2J). Writing the latter term in the standard form of a 
deformed potential (c.f. 19 ), the mean field part of H' reads as 

H'mf = H sp h - QJz - J^huJo/3 (q cos 7 - Q 2 sin 7) (9) 

where (3 and 7 are the deformation parameters. The self-consistent solution = 
shortly denoted as ) corresponds to an energy minimum of the energy surface 
E'(/3, 7) - (#')• 

We are aiming at the contributions of zero point quantum corrections to the 
moment of inertia stemming from the RPA vibrations about the above mean field 
solution. Accordingly, the quasi boson approximation (QBA) is applied to the 
Routhian H' , equation (jSJ), which is rewritten in the form 

H' = H' MF + ^ ^(Q m ) 2 -4e S (10) 

TCI 

Here, the third term is a separable residual interaction 



H ™= ^2^\Q m -{Q m )) (11) 



neglected in the mean field calculations. 
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Using the standard notation we refer to particle states (unoccupied single-particle 
orbitals) by subscript p and to hole states (below the Fermi level) by h. Defining the 
boson-like operators 



-*°P ( 12 ) 

the QBA means to treat the above operators l|12|) as exact bosons obeying the 
commutation relations 



b^b v 



0, 



(13) 



where the double index ph runs over all particle-hole pairs and is labelled by fi or v, 
respectively. In this approximation any single-particle operator F can be expressed as 

F= (F) +F {1) +F {2) (14) 

where the second and third terms are linear and bilinear order terms in the boson 
expansion, respectively. 

The final RPA Hamiltonian reads 

2 



^rpa - E E JbK - \ E - (&»)) 



(15) 



where = e' p — e' h is the energy of a particle-hole excitation with respect to the mean 

field part H^ F , equation @. We remind that in the QBA one includes all second 
order terms into the boson Hamiltonian such that (F— < F >) 2 = F^F^. 
The RPA Hamiltonian is diagonalised by solving the equation of motion 



^rpa ' d\ 



(16) 



for the phonon operators O x which are linear combinations of the basic bosons w„ and 
bp, equation (|12|l . The linear part of the Q-operators in terms of the boson operators 
b' and b has the following form 



Qr, 



lip, 



(17) 



Here q m ,^ are the single-particle matrix elements of Q m in equations I0HS|). 

The equation of motion fTSjl leads to the following determinant of the secular 
equations (c.f. [T^]) 



F(u x ) = det (R — — ) 
Ik 



where the matrix elements of R are given by 

Rkm{u\) = E ' 



involving the coefficients 



°p - 



E„ fc = 0,-l,2(l,-2) 



= 0,-1,2(1,-2) 



k = 0,-1,2(1,-2); rn = 1,-2(0,-1,2) 



The zeros of the function F 
F(u x ) = 0. 
yield the RPA eigenfrequencies ui\. 



(18) 
(19) 

(20) 
(21) 



Quadrupole correlations and inertial properties of rotating nuclei 



Since the mean field violates the rotational invariance, among the RPA 
eigenfrequencies there exist two spurious solutions. One "spurious" solution at u)\ = Q 
corresponds to a collective rotation, since [H', J±] = [H 1 , J x ±iJ y ] = ^fflJ±. The other 
solution with zero frequency is associated with the rotation around the z axes, since 
[H',J Z ] — 0. This spurious mode allows to determine the Thouless-Valatin moment 
of inertia, Jt.v. 1101113 EH, which can be calculated from the equations 

#BPA.*'*1 = ^ (22) 
J Jt.v. 

= i. (23) 



5 U Z 1 



Here the angle operator $ is the canonical partner of the angular momentum operator 

Jz- 

The total energy can be split in the mean field and the RPA contribution 

E ' = e 'mf + ^rpa + ~2 ( 24 ) 

with explicitly written the rotational "spurious" mode and the RPA contribution of 
non-spurious modes and of the spurious zero mode. The RPA contribution can be 
expressed as ^1 |H] 

b Rpa = 5 (X>*-E^ ■ (25) 

The exchange term ^3 [8] can be neglected which means to use the Hartree 
approximation. 

We recall that the mean field energy, the quasiparticle (particle-hole) excitations, 
the RPA eigenfrequencies, are calculated in the rotating frame. In other words, these 
quantities are functions of the rotational frequency f2 that is our external parameter. 
We remind that in the rotating frame the appropriate state variables are the Routhian 
energy E' and the rotational frequency £1 in contrast to the lab system where the 
appropriate state variables are the energy defined by the Legendre transformation 
E = E' +QI and the angular momentum /. While the operator relation dH' j d£l — —J z 
holds, the corresponding relation in the RPA order 

dE'/dQ = -(RPA\J Z \RPA) (26) 

is taken only as a reasonable approximation to the relation between the exact energy 
and the angular momentum. In fact, it is not only a very difficult task to calculate 
numerically the RPA expectation value (RPA\ J Z \RPA). This calculation sensitively 
depends also from the validity of the QB A itself which raises another problem ^5] . 
However, the use of the total Routhian and rotational frequency as the appropriate 
variables to determine the angular momentum of the system is consistent with our 
numerical calculations. 

The proper definition of the angular momentum within RPA is a known problem which 
was controversially discussed by several authors. We mention below some important 
points of these discussions. 

Assuming that the RPA expectation value of the angular momentum in the yrast state 
is 

(RPA\J Z \RPA) = I, (27) 
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Reinhardt proposed to use the equations l|26|) and (|27|> to define the quantisation 
condition for the angular momentum (see discussion in (221 where the rotation axis is 
chosen to be x). Since in the mean field approximation (J 2 ) = — dE' M¥ /d£l 00, we 
obtain from the above equations the following quantisation condition 

< J *>=' + 5 +f %* < 28 > 

According to Reinhardt, the smallness of the contribution of the RPA modes 
dE^ pA /dQ to the value of the angular momentum in the yrast state could be used as 
the validity criteria for the self-consistent cranking model. 

It should be pointed out that Marshalek raised the question about the validity of 
equation lt2*5|) (see discussion in J53j where the rotation axis is x as well). He assumes 
that the mean field value of the angular momentum in the yrast state is (J z ) = I 
which should be preserved in the RPA order too. He proposed to use the energy 
and the angular momentum as the appropriate variables to analyse the rotational 
properties. Considering the cranking Hamiltonian up to the second order of RPA, 
he proved that all conservation laws are restored (see Sec. 4.1 in ^2). Note that to 
prove all conservation laws it was important to keep the second order of the cranking 

(2) 

term ft J z in the cranking Hamiltonian as well. From this point of view, we would 
find the contribution of the second order term of J z operator to the expectation 
value (RPA\J Z \RPA) = (J z ) + {RP A\J^ ] \RP A) (see discussion in [El ED), which 
is inconsistent with Marshalek's assumption. Furthermore, the total energy in [T2*] 
and |2.'{| is a sum of the mean field energy Escc defined in the lab frame whereas 
the RPA correlation energy is defined in the rotating frame. Marshalek assumed that 
the RPA frequencies in the lab and rotating frame are the same. This leads to an 
inconsistency which becomes obvious if we consider a rotation about a symmetry 
axis. In this case the RPA states are characterised by the projection of the angular 
momentum upon the symmetry axis because [J 2 ,OjJ = \0\ with A being the value 
of the angular momentum carried by the phonons along the z axis. We thus obtain 

[#hpa, o{] = [H RPA - nJz, d\] = - xn)d{ = Lu x d{ (29) 

This equation explicitly demonstrates that the RPA eigenfrequencies in the lab frame 
(Da are different from the RPA eigenfrequencies in the rotating frame lo\ if the RPA 
modes carry angular momentum (see details and discussion in |25p. 

Thus, considering equations l|26l) and 1|27[) to be a suitable relation between the 
exact energy and the angular momentum and knowing the rotational dependence 
E' = E'(fl) of the total energy, it enables us to calculate the angular momentum 
/ = by the derivative 

/ = .7 ME + A Apa = -^-^p-i (30) 
and, subsequently, the dynamical moment of inertia as the second derivative 

= m < 3I > 

In practise the differentiations are numerically obtained with finite differences Af2 
as indicated in the above equations. Hereafter, we denote the angular momentum / 
calculated in accordance with equation l|3l)fl as Jrpa- 

The RPA solution u>\ — £1 that contributes to the total angular momentum by 
— 1/2, is very sensitive to numerical inaccuracies and to the limits of the configuration 
space in the calculation. In fact, at small but nonzero rotational frequencies the 
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spurious solution at zero and the solution at the rotational frequency get mixed up 
because of numerical inaccuracies. The rotational solution will then no longer be 
exactly at Q. Since the angular momentum is calculated as the derivative, these 
errors can lead to large erratic contributions to the angular momentum. We avoid 
such errors by explicitly finding the rotational solution and subtracting from the RPA 
energy. The missing term is then replaced by Ap. The rotational solution can be 
identified by its collective nature and large overlap with the J + ~ J x + iJ y operator 
to which it should be identical if the solution was found exactly. This method is 
applicable when the level density is low, which is satisfied when TiVl < 0.3 MeV. At 
larger Ml the rotational and the spurious solution decouple in a numerically stable 
way and the solution at the rotational frequency is found with a high accuracy. To 
calculate the RPA energy (|25|) it was crucial to apply the contour integral method 
developed in 

It is important to use a large enough configuration space in the calculations which 
include mixing of different A-shells. To reduce the number of TV-shells needed in the 
calculation we always transform the matrix elements into the stretched basis. 

3. The RPA quadrupole correlations 

In this section we focus our analysis on the contribution of the RPA correlations to the 
self-consistent mean field solution. Parts of these results where presented in |271 |5| . 

The self-consistent mean field Hamiltonian agrees with that of a rotating three- 
dimensionally deformed oscillator 



where U) x ,y,z are the three oscillator frequencies |T5) . 

We start with the version (i) of the QQ-model (A7V = 0, ±2 shell mixing is 
included) and without the volume conservation constraint. The self-consistent mean 
field solutions of the Routhian, equation JHJ, are found from the relations 



where, as usual, the contributions of exchange terms arc omitted. 

Two systems are considered: (a) the lighter nucleus N=Z=14 with the self- 
consistent deformation (e2 = 0.3,7 = 0, at O = 0) and (b) the nucleus N—Z=76 
defined with the corresponding set (e 2 = 0.3,7 = 0, at fi = 0) by choosing the 
appropriate strength k for each system. The largest TV-shell included in the calculation 
is N = 6 in (a) and N = 8 in (b). 

It will be seen that both systems behave qualitatively similar with respect to the 
rotation. By searching for the energy minimum for finite rotational frequencies the 
self-consistent mean field solutions are found. This minimisation makes also sure that 
the RPA will restore the rotational invariance. The mean field angular momentum 
for these systems calculated by means of equation (|3Uf) is shown at Figs. H El It is 




(32) 




(33) 
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-2 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

0.0 0.2 0.4 0.6 0.8 1.0 

hii [MeV] 

Figure 1. The angular momentum, in units of h, for a rotating deformed 
harmonic oscillator with (solid line) and without (dashed line) QQ (AN = 2) 
isoscalar RPA correlations as a function of Q. The self-consistent deformed 
oscillator is filled with 14 protons and 14 neutrons. 




5 i , , , , , , 

0.15 0.25 0.35 0.45 
hQ. [MeV] 

Figure 2. The angular momentum, in units of h, for a rotating deformed 
harmonic oscillator with (solid line) and without (dashed line) QQ (AN = 2) 
isoscalar RPA correlations as a function of Q. The self-consistent deformed 
oscillator is filled with 76 protons and 76 neutrons. 
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proportional to the rotational frequency f2. This means that our system even if the 
deformation is changing with the rotation gets an approximately constant moment of 
inertia |28| . 

Taking into account the RPA correlations, we obtain a substantial decrease of the 
total energy E 1 and also a dependence of the RPA angular momentum Jrpa on fi. It 
turns out that Jrpa oc f2. The dynamical moment of inertia is approximately constant. 
However, it is reduced (the reduction is 16% and 4% in N=Z=14 and 76, respectively) 
due to the RPA correlations, as seen in figure^ The RPA correlations cause relatively 
smaller effect on the moment of inertia in the heavier system than in the lighter one. 
This is because of the absolute ^-dependence of the -Erpa becomes slightly smaller 
and the mean field moment of inertia becomes much larger in the heavier system. 
The spherical oscillator part in our Hamiltonian, equation (J3J, can in these cases 
be replaced by any set of spherical single-particle energies. This does not affect the 
general properties of the RPA correlations. The result would not qualitatively change 
when including a spin-orbit splitting. 

Now the above systems are studied for the version (ii), i.e. by dropping the 
AN — ±2 matrix elements of the quadrupole operators Q m , equations Q-©- We 
recall that the quadrupole Hamiltonian with AJV = is commonly used for obtaining 
the self-consistent solutions which includes both non-zero spin-orbit terms as well as 
the QQ-forces [7]. This task needs much less efforts, since solutions can be found in 
a relatively small configuration space. Having in mind that the AN — 2 particle-hole 
excitations are energetically unfavored relative to the A7V = 0, one would expect a 
reliable results within the model. In addition, the effect of the higher lying AN = 2 
vibrations could be studied. The self-consistent mean field solutions of the Routhian 
in equation (JSJ are again found from the relations (K'>.'SII where this time the quadrupole 
operators are written in the stretched representation |2H1 EH ■ 

In this model the angular momentum is as above proportional to ft but the RPA 
correlation energy does not any more influence the moment of inertia as can be seen 
in figure 03 We conclude that an important part of the RPA correlations originate 
from the AN = 2 part of the QQ-forces. It is interesting to mention that the simple 
QQ-forces with AN = provide quite reasonable results with regard equilibrium 
deformations in the mean field calculations. 

An explicit calculation of the commutation relation 

#rpa.£]=0, (34) 

gives a stringent check of the RPA restoration of the rotational symmetry. In the 
numerical calculation of the N=Z=14 case above we found that equation (|34|l was 
fulfilled with an accuracy of 10~ 3 when the deformation parameters £2 and 7 were 
determined with the accuracy of 4 and 2 decimal figures, respectively. 

Next, we apply the volume conservation constraint and the double-stretched 
representation of quadrupole operators and compare the results with those obtained 
with the (normal) QQ-forces. We remind that the double-stretched representation 
is based on the condition such that the change in the density distribution must be 
accompanied with a change in the potential . 

For the triaxially deformed oscillator the self-consistency condition 

ul(x*)=u>t(yZ)=u J l(z>) (35) 
and the volume conservation constraint u> x uj y u) z — u>q yield a self-consistent residual 
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0.4 0.6 
hii [MeV] 



1.0 



Figure 3. The angular momentum, in units of h, for a rotating deformed 
harmonic oscillator with (solid line) and without (dashed line) QQ (A7V = 0) 
isoscalar RPA correlations as a function of Q. The self-consistent deformed 
oscillator is filled with 14 protons and 14 neutrons. 



quadrupole interaction |17l I18j 

H = H M p — H les 
where H' MF is defined by equation H32|) and 



^ , QmQr. 



(36) 
(37) 



m=0,±l,±2 

Here Q is the quadrupole operator written in doubly stretched coordinates 

Xi=Xi— , i = 1,2,3, (38) 

where u>i is one of the oscillator frequencies. The self-consistent k can be found as 

(39) 



_ 4tt MuA 

K = 



5 At 



From the self-consistency condition, equation l|35|l . follows that 
= 0, m = 0,±l,±2. 



(40) 



The Q operators do therefore not affect the mean field energy and can be used directly 
in the RPA formalism. 

In figure 01 results for the harmonic oscillator with 10 protons and 10 neutrons 
are presented. As for the cases studied above the QQ forces give a reduced angular 
momentum and lead to a delayed start of the rotational band at Tifl = 1.1 MeV (where 
the calculated angular momentum becomes positive) . The RPA correlations obtained 
with the use of the double-stretched residual interaction results in an increased 
dynamical moment of inertia in contrast to the results with the QQ-forces. The 
rotational dependence of the equilibrium deformations is presented in figure [S] The 
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0.0 0.5 1.0 1.5 2.0 
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4 
2 

-2 

0.0 0.5 1.0 1.5 2.0 
hft [MeV] 

Figure 4. The angular momentum, in units of h, for a rotating deformed 
harmonic oscillator with (solid line) and without (dashed line) isoscalar RPA 
correlations as a function of f2. Upper part: volume conservation, lower part: 
QQ-forces. The self-consistent deformed oscillator is filled with 10 protons and 
10 neutrons. Both systems have the same deformation at f! = 0. 

minimum calculated with the volume conservation constraint seems to be changing 
more slowly then with the QQ-forces. 

From figureE! which shows the mean field energy as a function of the deformation 
at different rotational frequencies with and without the volume conservation 
constraint, one can qualitatively understand the different results. The RPA 
contribution to the angular momentum is obtained as the negative derivative of the 
correlation energy. The QQ-forces give a very flat minimum that is changing relatively 
fast with the rotational frequency. This causes a large correlation energy that also 
is changing relatively fast with the rotational frequency. As a consequence, there is 
the large contribution from the RPA correlations to the angular momentum seen in 
figure A small contribution to the total angular momentum could be explained 
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Figure 5. The deformation parameters £2 and 7 of a rotating deformed harmonic 
oscillator calculated with the volume conservation constraint and with the QQ- 
forces. The self-consistent deformed oscillator is filled with 10 protons and 10 
neutrons. 
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Figure 6. Mean field energy as a function of deformation for the harmonic 
oscillator using the QQ-forces and the volume conservation constraint. For each 
method the curve is plotted for the non-rotating system and for ?if2 = 1 MeV. 



as the result of the increase of the rigidity of the potential with the increase of the 
rotational frequency f2. The flat potential energy surface is also the reason why we 
obtain a delayed start of the rotational band using the QQ-forces. With the volume 
conservation constraint one observes an opposite behavior. The minimum is much 
more distinct and is less affected by the rotation. In figure 0] a small increase in the 
dynamical moment of inertia can be seen which is consistent with the fact that the 
potential gets flatter with increasing f2 when using the volume conservation constraint. 
In figure [7| the mass dependence of the RPA contribution to the moment of 
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Figure 7. RPA contribution to the moment of inertia as a function of mass 
number A for N=Z harmonic oscillator with the volume conservation constraint. 
The plotted values are for £1 = (MeV). Upper panel: relative contribution. 
Lower panel: absolute contribution. 



inertia, as calculated with the volume conservation constraint and the double-stretched 
quadrupole residual interaction, is plotted. Even though the RPA contribution to 
is increasing with increasing a mass number A the relative contribution is decreasing. 
This is consistent with the simple picture that the mean field approximation works 
better in heavier systems. 

4. Summary 

The value of the angular momentum and the moment of inertia are generally dependent 
on the correlations induced by the shape vibrations. The size of the RPA correlations 
depends on the curvature of the mean field potential. Their influence on the angular 
momentum and on the moment of inertia is determined on how the potential curvature 
depends on the rotational frequency. This dependence is different in different mass 
regions and is determined by the degree of filling of the shell. The scale of the variation 
of the correlations is sensitive to the accuracy of the mean field approximation and to 
the size of the configuration space. 

The very flat potential of the QQ-forces leads to large effects of the quadrupole 
correlations upon the angular momentum and the moment of inertia. In fact, the large 
RPA contributions for the QQ-forces cause the delay of the rotation in contrast to the 
results obtained with the volume conservation constraint. In addition, the QQ-forces, 
without the volume conservation condition, may not produce any stable minimum 
above a certain deformation and it seems that the deformation in the considered cases 
is close to that limit. When comparing the results with and without the AN = 2 
mixing one notices that even when the two calculations provide similar mean field 
solutions, which is expected the RPA does yield significant differences. We 

conclude that these forces are not very realistic in describing RPA correlations in 
nuclear structure. 
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The RPA calculations with the volume conservation constraint and the double- 
stretched quadrupole residual interaction give a smaller contribution to the angular 
momentum. The RPA contributions to the moment of inertia is almost constant as 
a function of the rotational frequency but show a strong dependence on the mass 
number. Even though the absolute contribution is increasing with increasing the 
mass number, the relative contribution is decreasing from 11% at A=20 to less then 
2% at A=80. The RPA based on the double-stretched quadrupole residual interaction 
provides a more realistic result in contrast to the cases of the QQ-forces discussed 
above. 

Finally, we conclude that the volume conservation constraint is very important 
requirement when considering effective forces. Further analysis of the self-consistency 
between a mean field approximation and a treatment of the quadrupole shape 
vibrations for various realistic forces is needed to make reliable comparisons with 
experimental results. Another important question which has to be addressed in further 
investigations is how other multipoles such as monopole and octupole ones would affect 
the present results. 
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